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Motivation 


• TM Sensitive directions are 
not co-linear 

• Both TMs can be operated 
drag-free (within 
measurement band) 


• TM Sensitive directions are 
co-linear 

• Only oneTM can be 
operated drag-free, the 
other must be forced 


Potential noise source not 
present in LISA 
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Drift Mode 


Concept: 

• Concentrate TM2 forcing 
into small time windows, 
a.k.a. “kicks” 

• Allow TM2 to fly free 
between kicks 

• Compare measured TM 
noise with normal science 
mode to study force noise 
effects 


..Closed-! cop Kicks" 
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Data Analysis Challenge 
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• Periods of low noise (free-f lights) * 

separated by short periods of high 

noise (kicks) . a ; • 
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V 

• Must excise kicks 

L " 

• Care about frequencies down to 
I mHz but flights only last ~200s 

* i* 

• Need to combine information 

from multiple free flights to get 
low-frequency information ... 

• have to remove deterministic part 
of the free-flight 

• Desired Outputs 

• Residual spectrum or parameters 
that describe it 

• Parameters that describe 
Deterministic signal 
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Pre-processing 


• Read data from dynamical 
simulations into LTPDA 

• Remove any known cross 
couplings 

• Convert from displacement 
to acceleration 

• Identify & excise kicks 

• Remaining data stream has 
deterministic part plus quasi- 
periodic gaps. 
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Data Analysis Approaches 


• Apply window function to data that goes to zero during the gaps 

• Fill the gaps with randomly generated noise with the correct 
statistical properties 

• Compute and correct for spectral leakage effects caused by gaps 
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Gap-filling Approach 

Simultaneously do noise estimation and Bayesian parameter estimation of deterministic signals. 

N 

a(t) = n(t ) + 

i=l 

1 . Make a guess for underlying noise spectrum and coupling parameters 

2. Remove estimate of deterministic signal 

3. Fill gaps using “constrained Gaussian noise” 

4. Re-estimate spectrum 

5. Re-estimate deterministic pieces (Fisher matrix) 

6. Repeat 2-5 as necessary 

Does not require periodic gaps 
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Constrained Gaussian Noise 


• Want data with correct spectrum in 
frequency domain 

• In time-domain, this means correct 
2-point correlation function 
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— Input Target! 



Frequency [Hz] 


• Fill gaps with random data but adjust 
the mean of each data point based on 
correlation function and data outside 

gap 

• Caveat: In practice, need to truncate Cij 
at finite length...requires white or rising 
spectrum.Will need to pre-whiten. 
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Noise filling example 


Input Noise S„(f) = S 0 [1 + (///o)“] 





Frequency [Hz] 
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Noise filling example 


Add gaps (Is duration at 200s intervals) 




LPSD(yin) 
fitSpec[LPSD(yi 
LPSD(lteration I 
Best Fit 0 
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[m Hz 



Noise filling example 

First Iteration 



LPSD(yin) 
fitSpec[LPSD(yi 
LPSD(lteration ( 
Best Fit 0 
LPSD(lteration 
Best Fit 1 


Frequency [Hz] 
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Noise filling example 


Second Iteration 



Frequency [Hz] 
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f itSpec[LPSD(y i 

L PS D( Iteration i 

Best Fit 0 

LPSD(lteration * 

Best Fit 1 

— LPSD( Iteration ; 
Best Fit 2 
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Noise filling example 


Fourth Iteration 




LPSD(yin) 
fitSpec[LPSD(yi 
LPSD(lteration i 
Best Fit 0 
LPSD( iteration 
Best Fit 1 
LPSD( iteration : 
Best Fit 2 
L PS D( Iteration : 
Best Fit 3 
LPSD( Iteration < 
Best Fit 4 
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Convergence 

s„(f) = So [i + (///on 


Spectrum 

So 
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Time domain 


original 
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Combine the two steps - first try 
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Next Steps 


• Work on combined problem 

• Investigate other noise shapes, esp. red 
noise (pre-whitening?) 

• Apply to more realistic data (STOC?, 
SSM?) 
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